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ABSTRACT 

This paper proposes a model for post-seismic surface deformations by attributing them to 
lithospheric and asthenospheric viscoelasticity. The model predicts that the deformations due to 
lithospheric viscoelasticity depend on the ratio of the effective shear modulus acting long after 
the lithospheric viscoelastic relaxation to that acting immediately following the earthquake. 
While such deformations are generally smaller than those associated with asthenospheric visco- 
elasticity. they occur on a shorter time scale and may be in opposite direction to both the mo- 
tion occurring at the time of the earthquake and that occurring as the asthenospheric relaxation 
occurs. 
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POSTSHSMIC SURFACE Dl FORMATIONS 
DUE TO I ITHOSPHI RIC 
AND ASTHENOSPHI RIC VISCOELASTICITY 

There is considerable interest in the time dependent surface deformations following an 
earthquake. The description ol such deformations is important for understanding postseismic stress 
relaxation and redistribution, the occurrence of aftershocks, and in developing models of earth 
structure and rheology. Deformations that occur on a scale which is long compared to the time for 
elastic rebound in an earthquake are p esumably due to anelastie properties of the earth. We intend 
to model some of these anelastie properties by using linear viscoelastic elements to describe the 
earth’s rheology. An important principle associated with the linear theory is tnc correspondence 
principle which allows the time dependent viscoelastic deformation to be deduced by applying 
certain mathematical transformations to the solution for the associated elastic problem. The 
application of the correspondence principle to a homogeneous half-space reveals that when the 
viscoelastic properties in shear and dilatation are identical there is no additional surface deformation 
beyond the purely elastic displacement. Thus for time dependent postseismic displacement to 
occur with the context of a linear theory either there must be discontinuities in the properties of 
the half-space or the responses in shear and dilatation must be different. The latter view was 

adopted by Singh and Rosenman ( l l >74) who calculated surface deformations associated with a 
displacement discontinuity in a half-space which is elastic dilatational and Maxwell or Voigt 

deviatoric. Conversely. Nur and Mavko ( l l >74) considered a dislocation in an elastic layer lying 
over a viscoelastic half-space. Physically the elastic layer is associated with the earth’s litho- 
sphere and the viscoelastic half-space with the underlying asthenosphere. They considered both a 
thrust fault which was represented by an edge dislocation and a strike-slip fault which was repre- 
sented by a screw dislocation. Referring to the strike-slip case they found that the time-dependent 
displacement parallel to the direction of the initial fault motion is such as to increase the t.>tal 
displacement. The relative magnitude of the viscoelastic versus elastic disj lacement increases with 
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distance from the fault. Furthermore, the viscoelastic displacement increases as the depth of the 
screw dislocation approaches the lithosphere-asthenosphere boundary. The time constant associated 
with the motion is several years, this being the time constant for the viscoelastic response of the 
ast heliosphere. 

In addition to the viscoelasticity of the asthenosphere, there is evidence to suggest that a 
shorter time scale viscoelasticity can he associated with the lithosphere. Rock mechanics experi- 
ments, for example, suggest viscoelastic response times of I0 J - 10* seconds (Robertson, 1964). 

I pisodic creep events along the San Andreas fault have similar time constants (('rough and Burford, 
1977) and might be* associated with lithospheric viscoelasticity as might some aftershocks. 

With these thoughts in mind we propose the following model for postscismic deformations. 

The lithosphere's instantaneous elasticity and subsequent partial flow is 

irodeled by a standard viscoelastic solid with shear stress, o, -strain, e f relation 

(Fiugge ( 1967) 


dt) ° ■ ( < b + 1> ft) 
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The instantaneous shear modulus is The quantity r. = I ' is the retardation time at con- 

Pi ’ L 9o 

slant stress for the lithospheric layer. The asthenosphere is modeled as a Maxwell viscoelastic 
element with shear stress-strain relation 

y + Pr 


d \ de 

dt )° ' qj dt 
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Mere the shear modulus is — and the stress relaxation time for constant strain is r a = p. We 

di Pj q, 

shall have need ol the following parameters or = — — and 0 = t— - I . Numerically a has a 

V\ 4j "o"i 

.value near unity with the values 0.25 4.0 oeing used by Rybicki 1 1971); for 0 a value 0 = 3 is 

typical (I)ictcrich. 1972) l or the model ol the fault we choose the screw dislocation represen- 
tation of a strike-slip fault used by Rybicki ( 1971 ) and Nur and Mavko ( 1974). We shall find it 
convenient to assume r a > r, although this simplification is for physical clarity rather than 
mathematical necessity. 


The starting point for the analysis is the equation for the displacements associated witi. a 
dislocation. AU, in an el.«»;i«. layer with shear modulus overlaying a half-space with shear 
modulus /i 2 . With the coordinate system shown in Figure I . the only non-vanishing displacement is 
Uj where (Rybicki, 1971) 


U,(X, -0;X, >0)-f j- f*™ ' & f (frr£)""- [ 4 nV-Xv 


(3) 


To obtain the viscoelastic solution we replace fj, and f, i 2 by the corresponding Laplace transformed 

viscoelastic parameters jj, (s) and /T 2 (s) where s is the transform variable. We then replace All by 
AU 

since the solution is valid only for t > 0 and the Laplace transform lor this step function in time 
is j . The resulting modified equation I gives the Laplace transform for the viscoelastic deforma- 
tion. Taking the inverse transformation yields the desired time dependent displacements. To carry 
out this procedure we use the standard results (Flugge, 1967) 
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Then 
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The simple forms forS| and s 2 result from the assumption that r a ► r L . Taking the inverse Laplace 
transformation by the method of residues we find 
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where 


and 
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We are most interested in determining the surface displacements at three times. At time t = 0 we 
recover the instantaneous elastic displacements, at a later time r a ► t ► ft the deformation asso- 
ciated with viscoelasticity of the lithosphere has occurred and at yet a later time t > r. the 
deformation associated with viscoelasticity of the asthenosphere has also taken place. In these 
three limits equation 5 becomes particularly easy to evaluate and 


AU ) n 

Ujtt * 0) = ■ j 2 + tan 
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AU I n X, « 2X.I) I 
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a L I * ' 4n 2 H 2 -l) 2 + X, 2 \ 


AU| r Xi « _ 2X| D j 


(6) 


Ujft > T ) = — j- T+ tan* 1 -TT - 2 tan 1 ; } 

d I ~ 1 ■ 4n* H*-D 2 + X, 2 ( 


Notice that the deformations associated with the lithosphere viscoelasticity depend critically on the 
parameter B which in turn is strongly influenced by the quantity I + 0. This latter quantity is the 
ratio of the effective shear modulus for t > t l to that for t = 0. As this ratio increases so do the 
deformations. 
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As an illustration of the deformations predicted by these equations we consider the case a - 1, 
fi * 3, D * H/2. The results are presented in Figure 2. The displacements that occur on the time scale 
of the lithospheric viscoelasticity act to reduce the initial elastic displacements while those occurring 
on the time scale of the asthenospheric viscoelasticity increase the displacements. At first it might 
seem unlikely that the lithospheric viscoelasticity acts to reduce the initial elastic displacements. In 
fact, in an earlier paper concerned with the role of viscoelasticity in fault dynamics we argued that 
viscoelasticity could be responsible for forward post-seismic creep and aftershocks along a fault 
(Cohen, 1 ‘>78). The present behavior can he explained hy a simple one dimensional argument using 
the standard linear solid shown in Figure 3. Prior to the earthquake shear stress is slowly loaded on 
the system with the result that just prior to the earthquake the shear force is 

F: L 

Ma + *b 


where L is the initial stretch of the element. During the earthquake there is a sudden drop in the 
shear stress and stretch. Only the spring with elastic constant p a can respond on the time scale of 
an earthquake so immediately after the earthquake the shear force is 


F f=^a 


1 "W^h 


AX 


where AX is the displacement occurring in the earthquake. In the fault motion model there is 
sliding against friction during the earthquake and after the earthquake friction prevents further 
motion until the stress partially recovers to an amount exceeding the frictional resistance. In the 
present model for the motion off the fault there is no friction and the initial dislocation is relaxed 
by the sliding so that L = AX. Thus the elastic shear stress is negative immediately after the 
earthquake and viscoelastic flow acts to reduce the initial elastic displacement. We are, therefore, 
ignoring any additional forward displacement induced by post-seismic creep along the fault. 

The increased displacement occurring when the asthenosphere adjusts to the earthquake is the 
more usual viscoelastic rebound discussed by Nur and Mavko ( 1974 ) and is due to asthenospheric 
flow and associated drag on the lithosphere. 
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Although we have not yet investigated whether present geodetic or seismic data permits a 
definitive test lor the validity of the concepts presented here, there are a number ol potentially 
significant implications of the theory if it is correct. One of the most important is that surveys 
taken during the days to months following an earthquake may be misleading il interpreted in terms 
of the elastic rebound occurring at the time of the earthquake. Furthermore fa lure to achieve 
closure in a network survey might be attributed to the continued motion occurring as the network 
is traversed. On the positive side, repeated simultaneous determinations of the location of survey 
sites in the weeks following a major earthquake might produce much information about the earth's 
rheology and tectonic processes. 

Hie present theory has suggested that viscoelasticity in both the asthenospherc and litho- 
sphere may contribute to post-seismic surface deformations. A detailed quantitative model 
will require additional consideration of the three dimensional nature of the problem with the 
possibility that the dilational and deviatonc responses may be different. The model can be 
developed by applying the procedure used in this paper to the elastic deformation equations of 
Kundlc and Jackson 1 1977). The labor involved in such catenations is more arduous than that 
involved in those we have reported here; we are in the process of working out the details 
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FIGURE CAPTIONS 


Figure I. Model for Lithospheric and Asthenospheric Viscoelasticity. A Strike-Slip Fault is 
Modeled as a Screw Dislocation in the Lithospheric Layer. 

Figure 2. Computed Surface Deformations as a Function of Distance from the Fault for Various 
Times Following an Earthquake. 

Figure 3. Rheological Element for Standard Viscoelastic Solid. 



Figure I. Model lor Lithospheric and Asthenospheric Viscoelasticity. A Strike-Slip Fault 
is Modeled as a Screw Dislocation in the Lithospheric Layer. 
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Figure 3. Rheological I -lenient lor Standard Viscoelastic Solid 
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